Ultraviolet Probing of Quantum Crossbars 
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Ultraviolet (UV) scattering on quantum crossbars (QCB) is an effective tool for probing QCB 
spectral properties, leading to excitation of QCB plasmon(s). Experimentally, such a process corre- 
sponds to sharp peaks in the frequency dependence of the differential scattering cross section. The 
peak frequency strongly depends on the direction of the scattered light. As a result, ID — > 2D 
crossover can be observed in the scattering spectrum. It manifests itself as a splitting of single 
lines into multiplets (mostly doublets). The splitting magnitude increases with interaction in QCB 
crosses, while the peak amplitudes decrease with electron-electron interaction within a QCB con- 
stituent. 

PACS numbers: 78.67.-n, 78.30.-j, 73.90.+f 



I. INTRODUCTION 

Quantum crossbars (QCB) is a novel artificial nano- 
object which represents actually a double 2D grid 
formed by two superimposed crossing arrays of parallel 
conducting quantum wires 1,2,3,4 , molecular chains 5 or 
metallic single wall carbon nanotubes (SWCNT) 4,6,7,8 . 
Similar structures with the same crossbar geometry 
also arise naturally as, e.g., crossed striped phases of 
doped transition metal oxides^. The QCB mechanical 
flexibility together with the possibility of exciting some 
of its constituents (nanotube, a single wire) by external 
electric field and the existence of bistable conformations 
of others (molecular chamiS) makes QCB one of the 
most attractive architectures for designing molecular- 
electronic circuits for computational application 5 ^*!!. 
One more attractive property of QCB, its optical activity 
in a wide frequency region from infrared (IR) to ul- 
traviolet (UV) one, will be discussed in the present paper. 

From a geometrical point of view, QCB is, in a sense, 
an object with intermediate dimensionality. The two 
dimensional character of QCB as a whole, together 
with the one dimensional character of its constituents 
enables studying the crossover from a ID Luttinger 
liquid (LL) behavior to a 2D Fermi liquid (FL) behavior. 
Indeed, a single wire or nanotube possesses the LL-like 
spectrum 1 ^! 3 .. An array of parallel wires is still a LL-likc 
system (qualified as a sliding phased) provided the pure 
electrostatic interaction between adjacent wires is taken 
into account. However, inter- wire tunneling turns the 
electronic spectrum of an array to that of 2D FLi 4 *! 5 .. 
A system of two crossing arrays (QCB) manifests 
the same behavior. Being coupled only by capacitive 
interaction in the crosses, they have similar low-energy, 
long-wave properties characterized as a crossed sliding 
LL phas o 7 i 16 . Inter-array electron tunneling destroys the 
LL behavior at low energies and results in dimensional 
crossover to 2D FLii 7 ^. Here the control parameter 
which rules the crossover is the tunneling intensity. 

Outside the low energy region, the intermediate 
character of the QCB dimensionality leads to nontrivial 
spectral properties which cannot be treated in terms of 



purely II? or purely 2D electron liquid theory. QCB with 
only electrostatic interaction in the crosses possesses 
the LL zero energy fixed point. Far from this point, 
the system conserves a Bose character in the absence 
of interwire tunneling but the dimensionality of its 
Bose excitations (plasmon modes) becomes intermediate 
between ID and 2D in spite of the fact that they live in 
the 2D Brillouin zone (BZ) 1 ^. These QCB plasmons 
can be treated as a set of dipoles within the QCB 
constituents. In a single wire, the density of their dipolc 
momenta is proportional to the LL boson field 9(x) (x 
is the coordinate along the wire). In QCB, there are 
two sets of coupled dipoles. They form a unique system 
which manifests either ID or 2D properties depending 
on details of the relevant experiment. This results in 
a new dimensional crossover. It differs from the LL to 
FL transition mentioned above and is ruled by other 
control parameters (quasimomentum, energy, frequncy 
of an external field and so on). Some possibilities of 
observing such crossover in transport measurements 
(which give information about the nearest vicinity of the 
LL fixed point at (q,u),T) — > 0) were discussed in Ref0 
Other crossover effects such as appearance of non-zero 
transverse space correlation functions and periodic 
energy transfer between arrays ("Rabi oscillations") 
were studied in Ref^J Observation of these effects 
probes the QCB spectral properties outside the LL fixed 
point (i.e. well beyond the crossing sliding phase region). 

A rather pronounced manifestation of thess kinds 
of dimensional crossover is related to QCB response 
to an external ac electromagnetic field. The two main 
parameters characterizing the plasmon spectrum in QCB 
are the velocity v of plasmons in a single wire and the 
QCB period a (we assume that periods in the two basic 
directions are equal). These parameters define both the 
typical plasmon wave numbers Qq — 2ir /a and their typ- 
ical frequencies loq = vQq. Choosing v ss 8 • 10 7 cm/sec 
and a w 20 nm, according to Refs.|all3l|. one finds that 
plasmon frequencies lie in the far infrared (IR) region 
u> ~ 10 14 sec -1 , while characteristic order of plasmon 
wave vectors is Qq ~ 10 6 cto -1 . 



A very natural method to study the QCB spectrum 
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is by means of optical spectroscopy. The characteristic 
values of QCB frequencies and wave vectors determine 
two possible directions of such an experimental obser- 
vation. The first one is IR spectroscopy of QCB where 
the frequency of an external ac field lies at the same 
region as the QCB frequency. The second one is an UV 
scattering on QCB where the wave vector of a scattered 
field lies in the same region as that of the QCB wave 
vectors. 

The effectiven ess of IR spectroscopy of QCB were 
studied in Refs. I2fl2ll Here, the IR light wave vector 
k is three orders of magnitude smaller than the charac- 
teristic plasmon wave vector Q (the plasmon velocity 
v is much smaller than the light velocity c). Therefore, 
an infrared radiation incident on an isolated array, can 
excite plasmon only with lj — 0, or simply speaking, 
cannot excite plasmons at all. However in the QCB 
geometry, each array serves as a diffraction lattice for its 
partner, giving rise to Umklapp processes with Umklapp 
vectors say (riQo,0) with an integer n. As a result, 
excitation of plasmons in the BZ center with frequencies 
lj = nvQo is possible. 

To excite QCB plasmons with non zero wave vectors, 
an additional diffraction lattice with period A > a copla- 
nar with the QCB can be used. Here the diffraction field 
contains space harmonics with wave vectors 2-kM/A (M 
integer), that enables one to eliminate the wave vector 
mismatch and to scan plasmon spectrum within the BZ. 
In the general case, one can observe single absorption 
lines forming two equidistant sequences. However, in 
case where the wave vector of the diffraction field is 
oriented along some resonance directions, additional 
absorption lines appear. As a result, an equidistant 
sequences of split doublets can be observed in the main 
resonance direction (BZ diagonal). This is the central 
concept of dimensional crossover mentioned above with 
direction serving as a control parameter. In higher res- 
onance directions, absorption lines form an alternating 
series of singlets and split doublets demonstrating new 
type of dimensional crossover. The latter occurs in a 
given direction with frequency as a control parameter. 

One more version of IR QCB spectroscopy is related 
to the study of QCB placed onto a semiconductor 
substrateSSiSi. It occurs that a capacitive contact 
between the QCB and the substrate does not destroy 
the LL character of the long wave QCB excitations. 
However, the dielectric losses of a substrate are drasti- 
cally modified due to diffraction processes on the QCB 
superlattice. QCB - substrate interaction results in the 
appearance of additional Landau damping regions of the 
substrate plasmons. Their existence, form, and density 
of losses are strongly sensitive to variation of the QCB 
lattice constant. 

The IR based methods mentioned above are not very 
convenient from two points of view. First, as it was 
mentioned, one needs an additional diffraction lattice 



to tune the light wave vector and that of the QCB 
plasmon. Second, they probe QCB spectrum only in 
some discrete points. The alternative method of study- 
ing QCB spectrum by ultraviolet (UV) light scattering 
is the subject of the present paper. The advantages 
of this method are evident. It does not require any 
additional diffraction lattice. It probes QCB spectrum 
in a continuous region of wave vectors. Finally, its 
selection rules differ from those for IR absorption. This 
gives rise to the observation of additional spectral lines 
not visible in IR experiments. 

In this paper we formulate the principles of UV spec- 
troscopy for QCB and study the main characteristics of 
scattering spectra. The paper is organized as follows. In 
subsection III Al we briefly describe double square QCB 
and quantum numbers associated with its excitations. 
In the next subsection, III Bl we discuss light scattering 
on QCB and and present basic equations describing this 
process. The main results of the paper are contained 
in Section IIIII where we classify the basic types of the 
scattering indicatrices (angular diagrams of differential 
cross section) corresponding to various detector orien- 
tations. The results obtained are summarized in the 
Conclusions. Technical details are concentrated in two 
Appendices. Representation of QCB plasmons in terms 
of array plasmons is explained in Appendix^ Appendix 
iBlis devoted to the derivation of an effective QCB-light 
interaction. The basic formulas for differential cross sec- 
tion of light scattering are concentrated in this Appendix. 



II. LIGHT SCATTERING ON QCB 

A. QCB: Geometry, Hamiltonian, Spectrum 

A square QCB is formed by two periodically crossed 
perpendicular arrays of ID quantum wires or carbon 
nanotubes. The two arrays are lying on two parallel 
planes separated by an inter-plane distance d (see Fig. 
0. They are labeled by indices j = 1,2 and wires within 
the first (second) array are labeled by an integer index 
ri2 {n\). The coordinate system in Fig. ^ is chosen 
such that (i) the axes Xj and the corresponding basic 
unit vectors e, are oriented along the j-th array; (ii) 
the X3 axis is perpendicular to the QCB plane; (iii) the 
X3 coordinate is zero for the first array and d for the 
second one. A single wire is characterized by its radius 
ro, length L, and LL interaction parameter g. The 
minimal nanotube radius is ro ~ 0.35 nm2£, maximal 
nanotube length is L « 1 mm, and the LL parameter 
is estimated as g ~ 0.3^. Typical QCB period and 
inter-array distance are a w 20 run, d « 2 nm& so that 
usually ro < d < a < I. 

The QCB Hamiltonian is 

Hqcb = H 1 + H2 + H12. (I) 
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The first term H\ describes LL in the first array 

L/2 

Hi = -Q-yi / dx 1 lgn1(x 1 ,n 2 a,0) 

" 2 -L/2 1 
1 2 ) 

+ - (d Xl 0i (an, 7120,0)) ^, 
5 J 



(2) 



where (9i(x, y, z) and 7Ti(x, y, z) are two the conventional 
canonically conjugate boson fields of the first array. Its 
counterpart H 2 describes the second array and is ob- 
tained from Hi by permutation 1 <-> 2 and a replacement 
— > d in the arguments of the fields. 




FIG. 1: The scattering process geometry (all notations are 
explained in the text). 



The third term in Eq.QJ H12 describes the capacitive 
inter-array interaction. It has a separable form 



H12 



f L/2 



dxidx2 



L/2 

^ / xi - nxa \ / n 2 a - x 2 
x d xi 6x(xx,n 2 a, -d)d X2 6 2 (nia,x 2 ,0), 



(3) 



and is localized around the QCB crosses^. This form 
of inter-array interaction is actually a limiting form of 
an exact inter-array Coulomb interaction for ro/d <C 1 
taking into account electron screening in a single 
nanotube^. 

The Hamiltonian Q , J5J , © describes QCB in the 
absence of tunneling. In real QCB with typical param- 
eters mentioned above, an interwire tunneling is rather 



smali&. Coulomb blockade accompamying charging of 
adjacent tubes in a process of electron transition from 
one tube to another, suppresses these processes. Only 
electron-hole cotinncling without recharging is possible, 
and this cotunncling slightly renormalizes the coefficient 
Vq in Eq. ©. The corresponding relative correction is 
of order of 10~ 4 . 

The coordinates d and along the third spatial 
direction X3 enter the QCB dynamics as parameters. 
It is natural to expect that QCB possesses of 2D 
spectral properties which can be classified by 2D quasi- 
momentum. The basic vectors of the 2D reciprocal 
superlattice are Qo e j, j = 1,2, so that any vector of 
the reciprocal superlattice m is a sum m = mi + 1112, 
where xa.j = rrijQoej with integer rrij. The first BZ of 
QCB is \Qj\ < Qo/2 (see Fig. 0). An arbitrary vector 
Q = (QiiQa) of reciprocal space can be written as 
Q = q + m, where q = (<?i,<fe) belongs to the first BZ. 
QCB eigenstates are classified by quasi-momentum q 
and 2D band number S. 

However, the specific QCB geometry makes its spectral 
properties rather unusual. Consider an isolated array 1. 
Within the (xi,x 2 ) plane, its excitations are described 
by a pair of 2D coordinates (xi,n 2 a), i.e. a continuous 
longitudinal coordinate xi parallel to the ei direction, 
and its discrete transverse partner n 2 a parallel to e 2 . As 
a result, the longitudinal component Qi = q± + miQo 
of the excitation momentum changes on the entire axis 
—00 < Qi < 00 while its transverse momentum q 2 is 
restricted to the region \q 2 \ < Qo/2. Thus an eigenstate 
(plasmon) of the first array is characterized by the vector 
Qi = q + mi = (Qi,q 2 ) and the frequency 

Wi(Qi) =«|Qi| 

which depends only on the longitudinal component Qi 
of the momentum Qi. Similar description for the second 
array is obtained by replacing 1 <-> 2. 

In the case of "empty lattice" (QCB with infinitesimal 
interaction between arrays, Vq — > 0) the complementary 
second array produces infinitely small periodic potential 
which affects the first array. The vector mi becomes a 
reciprocal ID lattice vector and q becomes the 2D quasi- 
momentum. The first array excitations are labeled by 
this 2D quasi-momentum q and ID band number 



si 



2\Qi 



1 



(square brackets denote an integer part of the number). 
The corresponding dispersion law has the form 



w si (q) = wi(Qi 









. 2 . 



«|Qil 
+ (- 



; Qo 



As in the previous paragraph, by replacing 1 «-> 2, we ob- 
tain formulas describing eigenstates of the second array 
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Qi/Qo 




0.5 1 Qi/fio 

FIG. 2: Part of the reciprocal space. The small square 
Qj/Qo < 0.5 is the quarter of the first QCB BZ. High sym- 
metry lines are parallel to the coordinate axes. Resonant lines 
are parallel to BZ diagonals. The arcs coming over the point 
Bi, B2, -B3, B4, B5, correspond to different wave numbers of 
the excited plasmons: \Q/Qo\ = 0.3; 0.5; 0.6; 0.7; 0.8. 



of higher order of magnitude with respect to the small 
inter-array interaction J3J . 

Such a simple description fails at the points Q lying 
on two specific groups of lines in reciprocal space (see 
Appendix El f° r details). Nevertheless, even at these 
points QCB plasmons are represented as a finite su- 
perposition of array plasmons with coefficients of order 
unity. As the point Q moves away from a specific line, 
only one of these coefficients survives while all other 
decrease rapidly. So more complicated structure of QCB 
plasmons at the specific lines, presented in Appendix 
lAl turns to the standard description (jjjj within narrow 
transition strips which adjoin to the specific lines. Out 
of these strips, i.e. in major part of reciprocal space 
(including the nearest vicinities of the specific lines), 
a point Q completely determines the structure of all 
possible QCB plasmons with momentum Q. 



B. Light scattering 



embedded into an " empty lattice" . Thus an arbitrary 
vector Q = q + mi + m 2 = (Qi,Q 2 ) of the reciprocal 
space generates two vectors Qi = q + mi = (Qijffe) 
and Q 2 = q + m 2 = (<?i,Q 2 ). The "empty QCB" 
eigenstates are labelled by one of these vectors and its 
array number. Each one of the eigenstates propagates 
along the corresponding array. 

The dimensionless strength of the inter-array interac- 
tion of QCB with typical values of the parameters is 
small, 



The simplest process contributing to Raman-like light 
scattering is an annihilation of an incident photon and 
creation of a scattered photon together with a QCB plas- 
mon (Fig. In terms of initial electrons, this is in fact 
a second order process. Since the energies of incident 
and scattered photons significantly exceed the electron 
excitation energy in the nanotube, one may consider the 
emission/absorption process as an instantaneous act (see 
Appendix[BJ) . This process may be treated as an inelastic 
photon scattering accompanied by emission or absorption 
of a plasmon (Fig. Eh). 



hva 



0.007. 



(4) 



As a result, QCB eigenstates approximately conserve the 
same quantum numbers as the (unperturbed) eigenstates 
of an "empty QCB" . These numbers are array numbers 
j = 1,2, and the corresponding vectors Qj. In the gen- 
eral case, there are two QCB eigenstates associated with 
each wave vector Q of reciprocal space. In zeroth approx- 
imation, these QCB eigenstates coincide with the "empty 
QCB" eigenstates 



|l,Qi) = a\ Qi 
|2,Q 2 ) 



2,Q 2 



(5) 



where is an array of plasmon creation operators and 
|0) is the bosonic vacuum. The corresponding eigcnfrc- 
quencies coincide with the unperturbed frequencies 



u j (Q j )=v\Q j \, j = 1,2, 



(6) 



in the same approximation. In higher orders of perturba- 
tion, the r.h.s. of Eq. JSJ) includes additional terms like 

a iQi+ml0)' wnere m is an arbitrary reciprocal lattice 
vector. However, the coefficients before these terms are 



Let kj (k/) and £1^ = cki = c|k$|, (fi/ = cfc/) be 
momentum and frequency of the incident (scattered) 
photon. For simplicity we restrict ourselves by the 
case of normal incidence k^ — (0,0, fcj). Direction of 
the scattered wave vector is characterized by a unit 
vector n(ip, $), where tp and $ are polar and azimuthal 
angles in the spherical coordinate system with polar 
axis directed along e3 direction. The momentum of 
the excited QCB plasmon is Q = — K, where K is the 
projection of the scattered photon momentum k/ onto 
the QCB plane. This process is displayed in Fig. If 
the QCB plasmon frequency w(Q) coincides with the 
photon frequency loss lo = flj — fi/, then a detector 
oriented in a scattered direction n will register a sharp 
well pronounced peak at the frequency loss ld — w(Q). 
The frequency loss is much smaller than the incident 
and scattered photon frequencies. Therefore in what 
follows we use the same notation k for both ki and kf 
where it is possible. Scanning the frequency loss u>, (or, 
equivalently, the modulus of the scattered wave vector 
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X k, 



surface r = fo, and the vector of Pauli matrices <x is 



X k. 



k/- 



e 7 = — e2 sin 7 + e3 cos 7. 

The light wavelength is much larger than the nanotube 
radius, so the vector potential can be taken at its axis, 
A(zi,0,0). 



(a) 



(b) 



FIG. 3: (a) Second order diagram describing light scattering 
on QCB. Solid lines correspond to fermions, whereas dashed 
lines are related to photons. The vertices are described by the 
interaction Hamiltonian J7J. (b) Effective photon-plasmon 
interaction. QCB excitation is denoted by wavy line. 



kf) for a fixed detector orientation n, two or more (up 
to six for a square QCB) such peaks can be observed. 
The number of peaks and their location depend strongly 
on the azimuthal angle $ in the QCB plane. Scanning 
this angle, one can change the number of observable 
peaks. This is yet another manifestation of dimensional 
crossover mentioned above. 

The scattering differential cross section is smooth 
function of w. However, in zeroth approximation, only 
the peaks survive (smooth background is very low 
because of weak inter- array inetraction). Therefore in 
a chosen approximation, the scattering cross section is 
characterised by a number of peaks, their positions and 
intensities. Calculation of these quantities as functions 
of the frequency loss to and azimuthal angle $ for various 
fixed values of the polar angle ip, is the main goal of our 
study. 

To perform quantitative analysis, we should derive an 
expression for the scattering cross section. Details of all 
calculations are given in Appendix \B\ Here we present 
only the main steps of the derivation. We start from the 
Hamiltonian 



ev F 

n n i = 

c 



dx\drf 
2tt 



'Ipa' 



(7) 



which describes interaction between a single nanotube 
oriented e.g. along ei (the first array nanotube) and an 
external electromagnetic field. The field is described by 
its vector potential in the Landau gauge 



A = Aiei + A 2 e 2 + A 3 e 3 . 



(8) 



The indices a — a, b enumerate sublattices in a honey- 
comb carbon sheet, (r, 7) are polar coordinates in the 
(x 2 ,x 3 ) plane, ip a (xi,j) = tp a (xi, r , 7) and ip^ are 
slowly varying electron field operators at the nanotube 



Such form of a nanotube-light interaction leads to the 
following expression for an effective QCB-light interac- 
tion Hamiltonian 



Hi, 



V2e^ 
4fc he 



dxi x 



Here 



xd Xl 9i(xi,n 2 a)A 2 (xi,n 2 a,0) H 
4fc hc\ c J ^ J 

711 

xd x . 2 9 2 (nia, x 2 )Al(nia, x 2 ,d). 



Aj(r) = A(r) + (V2-l) ^(r)^, j = 1,2 



(9) 



are two effective vector potentials affecting the two arrays 
and Aj (r) are two (out of three) carthesian components 
of the full vector potential JSJ). The main object of our 
interest is the scaled differential cross-section ct(uj, n) of 
the scattering defined by 



'is: ($'(?)'•<"'•■>** <10 » 



da = L 



where L 2 is the QCB area and do = smipdipdi}. A stan- 
dard procedure applied to the Hamiltonian 0, leads to 
the following form of the scaled cross section 

^> n ) = iEi™°>pi 2 ^-^)- (u) 
p 

Here lu = f2,; — Q,f is the frequency loss and 

K; 



p.j. 



K 



(12) 



is the interaction Hamiltonian reduced to the subspace 
of QCB states. Summation is performed over all one- 
plasmon states \P). The vector Kj is obtained from K/ 
in the same way as the vector Qj is obtained from Q, 
and Pj;\t,\ i {<Pi'ff)i j = 1,2, are polarization matrices. In 
the basis (||,-L) they are 



2 sin d i cos -d 

2i cos 1? cos ip sin i? cos (p. 



P 2 (<p,#) = P 1 ^-+^, ( p 



(13) 



Equations i|ll|) - (|13|l serve as a basis for the subsequent 
analysis. 
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III. SCATTERING CROSS SECTION 
A. Cross section: Basic types 

According to Eqs. (| 1 1 11 . (|12|l . in order to contribute 
to the cross-section (|llf) in a fixed detector orientation 
at n, an excited QCB plasmon \P) must contain at 
least one of two single-array states |1, — Kj) or |2, — K 2 ). 
Analysis shows that there are five basic types of excited 
QCB plasmons depending on the location of the point 
K. Here we describe these types of plasmons and the 
corresponding structure of the differential cross section 
of scattering. In this description, beside the polar angle 
ip of the scattered photon wave vector kj, we use, in a 
sense, a mixed representation. It is based on the excited 
plasmon momentum Q = — K and the azimuthal angle 
d of the transverse component of the stattered photon 
wave vector K. 

i. General case: The point Q lies away from 
both the high symmetry lines and the resonant 
lines. This case is illustrated by the point A in Fig. 
There are two QCB plasmons \P), Eqs. ©, with 
frequencies (0, which contribute to the scattering. The 
differential cross-section is a sum of two peaks centered 
at these frequencies. After averaging over initial and final 
polarizations, it has the form 

a(uj, n) = Fi(n)5(uj - wi) + F 2 (n)5(uj - lu 2 ), 

where the functions 

Fi(n) = F(tp, #) = | sin^cosz?! x 

1 sin 2 ip cos 2 i? H — cos 2 w , 

L 4 4 J 

F 2 (n)= F^tf+I), 

describe universal angle dependences of the peak 
amplitudes. The functions Fi l2 are related to the cor- 
responding array plasmons. Each one of them vanishes 
when the scattered photon is perpendicular to the 
corresponding array. However these functions describe 
the scaled cross section. The absolute value amplitude of 

9 I e 2 \ 2 /vf \ 4 

each peak has an additional factor — — — — I — I . 

47rfca \ he J \ c J 

Strong electron-electron interaction in a nanotube 

corresponds to small values of the Luttinger parameter 

g and therefore suppresses the scattering cross section. 

ii. Inter-band resonance in one of the arrays: 
the point Q lies on a high symmetry line of only 
one array. This case is illustrated by the points C 2 ^ 
(1-st array) and D 2 % (2-d array) in Fig. [21 Consider 
for example point C 2 where Q\ — Qo/2, Q 2 ^ nQ /2. 
Here three QCB plasmons contribute to the scattering. 
The first is a |2, Q 2 ) plasmon, (JjJJ, with frequency uj 2 , Eq. 
© . The other two are even or odd superpositions of the 
1-st array states (Eq. 1|A1J1 with j = 1) of the two first 
zones with eigenfrequencies l|A2(l . Due to weakness of the 
inter-array interaction, three peaks of the scattering cross 



section form a singlet oj 2 and doublet u>i g ,uji u . After 
averaging over initial and final polarizations, the cross 
section has the form 

a(u>,n) = -Fi(n)[8(u - uig) + 5(u - u>iu))] 
+ F 2 (n)5(v - uj 2 ). 

iii. Inter-band resonance in both arrays: The 
point Q is a crossing point of two high symmetry 
lines away from all resonant lines. This case is illus- 
trated by the points C4 and D4 in Fig. [2 Consider for 
example point C4. Here Qi = Qo/2, Q 2 = 2Q /2, and 
four QCB plasmons contribute to the scattering. The 
first pair consists of even and odd superpositions of the 
1-st array states of the first and the second bands. These 
states and their frequencies are described by Eqs. l!A5|) . 
(IA6I) with j = 1 . The second pair consists of the same su- 
perpositions of the 2-d array states from the second and 
third bands and is described by the same equations with 
j = 2. As a result, four peaks, which form two doublets 
(|A6I) . j = 1,2, can be observed. After averaging over 
initial and final polarizations the cross section is 

= -Fi(n)[5(uj - uj lg ) + 5{uj - uj 1u ))} 
+-F 2 (ii)[S(uj - uj 2g ) + S(uj - uj 2u ))}. 



iv. Inter-array resonance: The point Q lies only 
on one of the resonant lines away from the high 
symmetry lines. This case is illustrated by the points 
£?i_3,5, E, and F in Fig. [3 Here the QCB plasmons 
which contribute to the scattering are two even and two 
odd superpositions of the first and second array states 
(|A3I) whose eigenfrequencies form two doublets (|A4I) . As 
in the previous case, the scattering cross section contains 
four peaks which form two doublets. After averaging 
over initial and final polarizations the cross section is 

<r(w,n) = i J Fi(n)[<5(w-cj g (Q 1 )) + <5(w-w n (Qi))] + 

^F 2 (n)[S(cu - u g (Q 2 )) + S{u; - u u (Q 2 ))}. 

Thus, the inter-array splitting is proportional to the 
central small parameter of the theory, tf> (see Eq. Q). 
For the set of parameters described in the beginning 
of subsection III Al the interband splitting defined by 
Eq. i|A2(l . is five times smaller because it contains an 
additional factor tpa/rQ. 

v. Inter-array and inter-band resonance: The 
point Q lies at the intersection of two resonant 
lines. There is only one such point B4 in Fig. In 

the general case, where the parameter n / for both 
crossing resonant lines (the point B4 is not the case), the 
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QCB plasmons involved in Raman scattering form two 
quartets. The first quartet consists of four symmetrized 
combinations (|A7|) of the single-array states. QCB eigen- 
states for the second quartet are obtained from these 
equations by replacing 1 <-> 2. The corresponding eigen- 
frequencies are described by Eqs. IjASjl . (|A9(I . The scat- 
tering cross section in this case contains six peaks, and 
two of them are two-fold degenerate 

cr(w,n) = ^Fi(ii)[2S(uj - w g>u / 9 (Qi)) + 
S(oj - uj uu (Qi)) + S(uj - w„ ff (Qi)] + 

S(oj - uj uu {Q 2 )) + S(u - u ug (Q 2 )]. 

The point B4 in Fig. lies on the main resonance 
line with n = 0. Here |Qi| = IQ2I, the frequencies of 
both quartets coincide, and the scattering cross section 
contains one four-fold degenerate peak and a symmetric 
pair of its two- fold degenerate satellites. 

This classification of all types of excited plasmons 
enables us to describe completely the UV scattering on 
QCB. 



B. Scattering indicatrices 

In this subsection we describe the results of the 
scattering process with the help of a family of scattering 
indicatrices. 

To explain the indicatrix structure, we start with 
some preliminary arguments. Consider the case where 
the detector is tuned to the frequency Q f and is oriented 
in the direction n. These parameters uniquely determine 
a point K = Q. There are two ways of scanning QCB 
plasmons. The first way is to scan through the polar 
angle ■&. The corresponding points K in Fig. [21 form an 
arc with radius K — (flf /c) sin (p. The second way is 
to tune the detector frequency f2/. We are interested 
in frequency loss of order of the plasmon frequencies 
within the two - three lowest bands. This loss is 
of order of 0/(u/c) <C flf. Therefore in this case the 
point K remains in its place (with a very good accuracy). 

Each scattering indicatrix corresponds to a circular 
arc in Fig. and the structure of this indicatrix 
is completely determined by the arc radius K. The 
indicatrix represents a set of curves displayed in polar 
coordinates with polar angle, which coincides with the 
azimuthal angle & used above, and (dimensionless) 
radius u)/(vQq), where to is the frequency loss. Each 
point of the indicatrix corresponds to an excitation of 
a QCB plasmon and therefore to a sharp peak in the 
scattering cross section. The number of peaks depends 
on the polar angle. Scanning the azimuthal angle i9 
results in changing the number of peaks. This is one 



more example of dimensional crossover in QCB (see 
Refl2(l for similar effects in IR spectroscopy) . 

We start with the case of smallest radius 
k sin if = — 0.3Qo (arc AB\ in Fig. 0). Here all 
points beside the point B\ are points of general type 
(i). Each one of them, e.g. point A, corresponds 
to excitation of two plasmons in the two arrays and 
therefore leads to two separate peaks in the scattering 
spectrum (see Fig. 0J. The peaks corresponding to the 
point A in Fig. lie at the ray defined by the angle 
0A- The point B\ corresponds to inter-array resonance 
(iv) and in this direction a split doublet can be observed. 




FIG. 4: Positions of the scattering peaks for |_K"/Qo| = 0.3. 
The doublet in the resonance direction (point B\ in Fig. |5J 
is well pronounced. 

In the next case fcsint^ = — 0.5Qo ( arc C\B 2 D\ in 
Fig- EJi as in the previous one, in general directions 
one can observe two single lines which form a split 
doublet in the resonant direction B 2 . However at 
the final points C\ and D\ the arc touches the high 
symmetry lines. Here the low frequency single line 
vanishes (there is no scattering at Q = 0) while the 
high frequency line transforms into a doublet because of 
an inter-band resonance (ii) in one of the arrays (Fig. [SJ. 

Further increase of the arc radius ksinip = — O.6Q0 
(arc C 2 B$D 2 in Fig. [2Jl leads to appearance of two points 
D 2 and C 2 where the arc intersects with high symmetry 
lines (BZ boundaries). Each of these points generates an 
inter-band resonance doublet, which coexists with the 
low frequency single peak (Fig[SJ). 

In the case ksimp — — y/2Q /2 the corresponding arc 
includes the BZ corner B4. This is the point of a double 
inter-array and inter-band resonance (v). Moreover, here 
the two quartets described above coincide. Therefore, 
there are three lines in Fig. [7| The low-frequency line as 
its high-frequency partner is two-fold degenerate while 
the central line is four-fold degenerate. We emphasize 
that each quartet manifests itself in three lines, contrary 
to the IR absorption, where selection rules make two of 
them invisible. 
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FIG. 5: Positions of the scattering peaks for \K/Qo\ = 0.5. 
Two doublets appear at the BZ boundaries (points Ci and 
Di in Fig. 0. 




FIG. 6: Positions of the scattering peaks for \K/Qo\ = 0.6. 
Two doublets at the BZ boundaries (points C2 and D2) in 
Fig. |5J are shifted from the high symmetry directions. 



The last case fcsint,!? = — O.8Q0 demonstrates one more 
possibility related to inter-band resonance simultane- 
ously in two arrays (iii). Each point E and F generates 
(in the corresponding direction) two doublets describing 
the inter-band splitting in different arrays (see Fig. [SJ. 



IV. CONCLUSIONS 

In conclusion, we studied inelastic UV Raman scat- 
tering on QCB. We derived an effective Hamiltonian for 
QCB-light interaction which is expressed via the same 
bose variables that the QCB itself. With the help of 
this Bamiltonian we calculated differential scattering 
cross section as a function of detector orientation and 
scattered frequency. Scanning these parameters, one can 
observe a set of sharp peaks in the scattering spectrum. 
A number of peaks and their positions strongly depend 
on the direction of the scattered wave vector. This 
results in a dimensional crossover. It manifests itself in 
the splitting of the peak frequencies and therefore in 
appearance of multiplets (mostly doublets) instead of 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 wvQo 

FIG. 7: Positions of the scattering peaks for \K/Qo\ = V2/2. 
Resonance triplet corresponding to point B4 in Fig. (two of 
four frequencies remain degenerate in our approximation) . In 
IR experiments only one of the triplet components is visible. 




FIG. 8: Positions of the scattering peaks for \K/Qo\ = 0.8. 
Two pairs of doublets appear corresponding to excitation of 
two pairs of plasmons in the two arrays (points E and F in 
Fig.0. 



single lines in the scattering spectrum. 

The sizes of peak splitting are determined by the 
nature of interaction which lifts the corresponding degen- 
eracy. In the case of initial inter-array degeneracy, the 
splitting is proportional to the dimensionless interaction 
strength (J2J. The inter-band splitting is proportional 
to the square of this interaction. For a chosen QCB 
parameters it is less than the inter-array interaction 
strength in spite of an additional large multiplier a/rQ. 
In all cases the splitting increases with increasing the 
interaction in QCB crosses. The peak amplitudes are 
proportional to the Luttinger parameter g in a single 
nanotube. Therefore strong electron-electron interaction 
suppresses the peaks. 

The effectiveness of UV scattering is related to the 
possibility of changing continuously the excited plasmon 
frequency. Due to other selection rules, some lines invisi- 
ble in the IR absorption spectrum, become observable in 
the UV scattering. UV scattering spectroscopy enables 
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one to restore parameters describing both interaction in 
QCB crosses and electron-electron interaction in a single 
QCB consituents. 

Our st udies o f optical properties of QCB (this paper 
and Refs 20 21) show that these nanoobjects possess 
unique combination of optical spectra. Firstly, they are 
active in IR and UV frequency range. Secondly, they 
may be observed in various kinds of optical processes, 
namely direct and indirect absorption, diffraction, 
energy loss transmission, and Raman-like spectroscopy. 
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APPENDIX A: QCB PLASMONS: KINEMATICS 
OF INTERFERENCE. 

As it was mentioned in subsection III Al standard 
description (JjJJ, 0, of QCB plasmons in terms of 
array plasmons fails at points Q lying on two specific 
groups of lines in reciprocal space. Nevertheless, QCB 
plasmons even at these lines can be considered as finite 
combinations of array plasmons generated by the point 
Q by more complicated way. Below we consider var- 
ious types of the inrterference of array plasmons in QCB. 



are even or odd superpositions of the 1-st array states 
from the two first zones with eigenfrequencies 

^i ff (Q) = HQil, 

w lu (Q) = (l-f^jvlQrl (A2) 

The case j = 2 is described similarly after change 1 <-> 2 
through all this paragraph. 
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FIG. 9: QCB inverse space. The first BZ is bounded by a 
sooid square 



The first group of specific lines is formed by the high 
symmetry lines Qj = pjQo/2 with pj an integer (lines 
parallel to coordinate axes in Fig. The lines with 

Pj = ±1 include the BZ boundaries. In an "empty 
lattice" approximation, these lines are degeneracy lines 
which separate pj-th and (pj + l)-th bands of the j-th 
array. Inter-array interaction mixes the array states, lifts 
the (inter-band) degeneracy and splits corresponding 
frequencies. As a result, the QCB plasmons related 
to a point Q lying ona high symmetry line (e.g. all 
C, D points in Fig. are built from the array plas- 
mons associated not only with the points Q1.2 but 
also with the symmetric (with respect to coordinate 
axes of a reciprocal space) points Qy = (— Qx,Q2) an d 
Q2 = (91, -Q2) (see Fig. EJ. 

Consider for definiteness the case j ' = 1 and assume 
first that the ratio 2Q2/Q0 is non- integer. In this case, 
the point Q generates three QCB plasmons. The first one 
of them is a |2, Q2) plasmon (see Eq|SJl with frequency 
u)2- The two others 

|l, fl> Q0 - -J=Ql,Qi) + |l,Q T )), 

|l,«,Qi) = -4(|i,Qi)-|i,Qt» (Al) 



The second set is formed by the resonant lines defined 
by the equation Q\ + rQi = nQo, where r = ±1 and 
n are two integer parameters determining the line. In 
Fig. [21 resonant lines arc square diagonals and lines 
parallel to them. In an "empty lattice" approximation, 
these lines are also degeneracy lines. However this 
degeneracy has more complicated nature. Here, array 
plasmons associated not only with the point Q, but 
also array plasmons associated with the dual point 
Q = (— rQ%, —tQx) (see Fig. HOf) are involved in the res- 
onance. Inter-array interaction mixes degenerate modes, 
lifts the degeneracy, and splits degenerate frequencies. 
As a result, the QCB plasmons related to a point Q 
at the resonant line (e.g. all B,F,E points in Fig. 
are built from the array plasmons associated not only 
with the points Qi,2 but also with the symmetric (with 
respect to one of the bisector lines of the coordinate 
system of a reciprocal space) points Q : = (— rQi, —rqa) 
and Q 2 = {—rqi, —rQ-i)- Moreover, if the point Q lies at 
the intersection of two resonant lines (like the point B4 
in Fig, the QCB plasmons include also array states 
associated with the Qy points. 

In the general case, the point Q lies only on one of 
the resonant lines away from the high symmetry lines. 
Here the point Q generates two pairs of even and odd 
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superpositions of the first and the second array states 
1 



\g/u,Qi 



v/2 
1 



(|l,Qi)±|2,Q 1 )) 



\g/u,Q 2 ) = ^=(|2,Q 2 )±|1,Q 2 )) 
whose eigenfrequencies form two doublets 



w 9 /„(Qi) = (l±|J«|Qi| 
„iQ,.) ' (l±| ) <1<7->i 
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FIG. 10: Wave vectors involved in the resonance interaction. 
Resonant lines are displayed by solid lines. 



The second group consists of the crosses of two reso- 
nant lines (point B4 in Fig. [2J. These point are always 
the crosses of two high symmetry lines. Here the QCB 
plasmons, generated by the point Q, form two quartets. 
The first quartet is really generated by the point Qi. It 
consists of four symmetrized combinations of single-array 
states 



\g,g,Qi 



\g,u,Qi) = 77 



ks,Qi) 



\u,u, Qi 



(|i,Qi) + |i,Qr)- 

|2,Q 1 ,) + |2,Q T >); 
(|1,Q!> + |1,Q T > - 
|2,Q 1 )-|2,Q T >); 

Qi) + |1, Qr) - 
|2,Q 1 )-|2,Q T )); 

(|1,Q!> - |1,Q T > - 
|2,Q 1 ) + |2,Q T 



(A7) 



QCB eigenstates of the second quartet are generated by 
the point Q 2 . They are obtained from equations l|A7(l by 
replacing 1 «-> 2. Even array eigenstates are degenerate 
with frequencies 



w g, 9 /ti(Ql,2) — 

while the odd array eigenstates are split 

Vu,g/u(Ql,2) = (1 ± 0)^1,2- 



(A8) 



(A9) 



The crossing points of resonant lines can also be di- 
vided into two groups. The first group is formed only 
by the crossing points of two high symmetry lines (e.g. 
points -D4, C4 in Fig. |2J). Here each of the generated vec- 
tors and Q 2 corresponds to a pair of array eigenstates 
belonging to two adjacent ID bands. Thus an inter-band 
mixing is significant in both two arrays. The correspond- 
ing QCB plasmons are even or odd combinations of j-th 
array plasmons (j = 1,2) 



\j,g,Qj) = 
\j,u,Qj) = 
with eigenfrequencies 

Uj,g(Qj) 



b'.Qj) 
QoPj 



\j,Qj. 



(A5) 



T a\ Q p 3 

1 V — — i 



ft 

ro 



(A6) 



Equations IjAlfl - l|A9|l exhaust all cases of QCB 
plasmons, generated by a point Q 2 lying on a specific 
line, via the interference of array plasmons. 



APPENDIX B: DERIVATION OF T HE BA SIC 
EQUATIONS OF SUBSECTION ITTEl 

Light scattering on QCB is described by equations 
- <|13|) . In this Appendix we briefly explain the main 
steps which lead to such a description. 

1. Nanotube-light interaction, Eq. ©. We start 
with a consideration of a single nanotube of the first ar- 
ray, interacting with an external electromagnetic field. 
The characteristic time of all nanotube energies, includ- 
ing Coulomb interaction, is of order of an inverse plas- 
mon frequency. The scattering process occurs during 
much shorter time interval, which is of the order of an 
inverse photon frequency. Hence, Coulomb interaction is 
irrelevant for the scattering processes. This enables us 
to restrict ourselves to the kinetic part of the nanotube 



11 



Hamiltonian. Within the k — p approximation, this part 
is given by24 



Vf 



dxidj 
2tt 



ipl(xi,r ,j)Kk- cr a . a >ip a i(xi,r ,j). (Bl) 

In the presence of a magnetic field, one should add eA/c 
to the electron momentum operator frk. An additional 
part of the Hamiltonian <|Blll exactly coincides with the 
nanotube-light interaction Hamiltonian^. 
Hamiltonian (|B1|1 is diagonalized by a two step canonical 

a>r /v F 
2 



kr 



FIG. 11: The low-energy part of the band structure of a metal- 
lic nanotube. Lines 1,3 (2,4) correspond to excitations with 
p — +1 (p = —1). Lines 1,2 (3,4) correspond to excitations 
with orbital quantum number m = (m = ±1). 



transformation. The first one is Fourier transformation 

c a .k, m e lkxi+m \ (B2) 



k,m 



where the orbital moment m is restricted by the condi- 
tion \m\ < niQ = [7rro/ao] due to the finite number of 
honeycomb cells along the nanotube perimeter is 
the lattice constant). The second rotation is defined as 



~y= (Cakrn + pe^"' C bkm ) , p = ±, 



where 



cos< 



krn 



sin ( 



y/ikro) 2 + m 2 V ( kr o) 

As a result, Hamiltonian takes the form 

p,k,m 

with eigenfrcqucncies 



uj p (k 1 m) = pvf\ k 2 + 



(B3) 



(B4) 



2. QCB-light interaction, Eq. ©. Substituting into 
Eq. the scalar product A • <x in the form 



(A- a) 



A~ 
A+ 



where A^ — A\ ± iA^, we can write the nanotube-light 
interaction as 



H,, 



evF 
c 



dxidj 



i ^^ A A-^ B +h.c.y 



We are interested in an effective QCB-light interac- 
tion Hamiltonian obtained in second order of perturba- 
tion which describes transitions between initial states \i) 
and final ones |/). Initial states are one-photon states of 
the electromagnetic field and the electron ground state of 
the nanotube whereas final states consist of one photon 
and an electron above the Fermi level. It will be seen later 
that just these states form a one-plasmon array states. 
The energy of incident photon E = hck is much larger 
than the excitation energies of the nanotube. Therefore, 
absorption of the incident photon and radiation of the 
scattered photon occur practically without retardation. 
All these result in the interaction Hamiltonian 



hint — 



he V c / 



dxidj f dx^d'y' 



2tt 



2tt 



T 1 - 



(^O^YM-o^o.ojiv) x 

(v\ij b (x'i,7')^t(xi,j)\v) x 
(«|A+(a: 1 ,0,0)^ o (a;i,7) + 

a <-> b , 



(B5) 



which corresponds to the diagram shown in Fig. [3^ 
(there is no photons in an intermediate state). 

Consider now the matrix element 
(v\tp a (x' 1 ,Y)'4'l l (xi,-f)\v) which enters this Hamilto- 
nian. Due to our choice of initial and final states, only 
diagonal elements with respect to both virtual states 
\v) and sublattice indices a survive. In the c a ^,m 
representation, (|B2|) they have the form 

{v\rj} a {x'x,^)i>i{xx,i)\v) = Y x 

k,k f m,m' 

(«|ca,fe,™cl ife ,, m »e lfca:i - 4fc ' a; ' 1 +^-™'^. 

Internal matrix elements on the r.h.s. of the latter equa- 
tion are 

(v\c a , k , m cl k , m ,\v) = \8 k<m 8 k , <m , 

(here the symmetry property n{E_(k,m)) + 
n (E + (k,m)) = 1 is used). Therefore 

(w|V>a(a;i,7')V'a0ci,7)k) =7r<5(2:i -x' 1 )S'(7-7'), 
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where 



1 



mo 

5(7) = f- E e "" 7 - 

Z7T * — ' 

Tci— — mo 



Thus, the interaction (|B5|> takes the form 
- 2 ' "■■ ■ ' r d.r ihih ' 



_ f VF \ 
he \ c J 



Ank 



S(l - 7') x 



i-( I i,0,0)#( I1 ,0 1 0)^i(i 1 , 7 ')i(ii,7). 

a 

The field dependent factor here is 

A-(x 1 ,0,0)A + (x u 0,0) = Al(x x ,0,0) + 
A 2 (x!,0,0) sin(-f A - 7')sin(7A - 7), 



(B6) 



and (Ai, A,ja) are cylindrical components of the vec- 
tor potential A (JSJ). Taking into account the angular 
dependence of the field (|B6|) . we can omit it in the elec- 
tron operators. Indeed, according to Eq. I|B4(1 (see also 
Fig. Ill|) in an energy-momentum region where we work, 
the m = 1 spectral band with nonzero orbital moment 
is separated from that with m — by an energy of or- 
der hvp/fQ, which is much larger than the QCB plasmon 
energy. Keeping only the zero moment field operators 
which form the electron density operator 

E V4( x i)V>a(zi) = P( x i) = V2d xi 9(xi). 

OL 

and integrating over 7,7', we obtain an interaction 
Hamiltonian in the form 



V2< 



( VF 



4fc he V 



' J d Xl d xi e{ Xl )Al{x u Q,0), (B7) 



where 



Ai = A + (V2 - l)Aiei. 



Straightforward generalization of this expression to the 
QCB case leads exactly to the Hamiltonian IjjJJ. 



3. Polarization matrix, Eq. (|13|h To study the 
scattering process, we should modify the last expression 
for the interaction Hamiltonian. To proceed further, we 
define Fourier transforms Oj,Q j of the bosonic fields 



i{n 1 a,x 2 ) 



Hi 

—iQxXt —iq 2 n 2 a 

=E^- 

ATT t-^l ' 



xe 



Q 2 x 

Q2 

iqin 1 a— iQ 2 x 2 



(B8) 



Here N = L/a\s the number of QCB cells in both direc- 
tions. The electromagnetic field also can be expanded in 



a sum of harmonics with a wave vector k and polarization 
A = |U, 



A(r) = E n kA-4kAe 



ikr 



(B9) 



kA 



The polarization vectors 

ik 2 ei ik\e 2 



n k || 



n k j 



VW+M 



VW+M 



fcifc 3 ei k 2 k 3 ei 



2 -L- k 2 
ft, 2 



k \jk\ 7 k<2 k \jk\ ~\~ k^ 



e 3 , 



are normalized, |nkA| — 1 5 and satisfy the orthogonality 
conditions, n k A ■ k = n k || • n k j_ = 0. The field oper- 
ators AkA satisfy the condition A^ x = A-^\, so that 
At(r) = A(r). 



Substituting equations l|B8|) . and (|B9fl into the Hamil- 
tonian @, we obtain 



Hint 



.\f2NL e 2 /v F \ 



4k he 



(?) S £ 



k,k',Q j,A,A' 



k' k 

A 7 ' fc 



Here 



Pi: 



k' k 



is the polarization matrix, Q = q + m, and 
«j,k,A = n k A + [V2 - l) (nkA • 



(BIO) 



(Bll) 



In the case of normal incidence, Eqs. ijBllfl . (|B12|) result 
in the form (|13() of polarization matrix. 

4. Scattering Hamiltonian, Eq. (1121) . As a next 
step, we express the Fourier transforms of the Bose fields 
9 via creation (a^) and annihilation (a) operators of the 
array plasmons 



The electromagnetic field amplitudes A kj A should also be 
expressed via photon creation (c^) and annihilation (c) 
operators 



A k , A (i) 



he 
2Vk 



Substituting these expansions into Eq. l|B10fl we obtain 
the final form of the effective interaction. In the case of 
normal incidence (K = 0), this interaction is written as 



Hi, 



V \2ck~) X 

E E *v,aCko4, iV *,. 



(B12) 



k,k'i,A,A' 
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where h\>,\(K') is the Hamiltonian (I12fl . where A/, A;,K 
are replaced by A', A, K'. 

5. Scattering cross section, Eqs. (Il()[) - (II 111 . Stan- 
dard procedure based on the Fermi golden rule leads to 
the following expression of the differential scattering cross 
section per unit QCB square 



1 fVk 



1 da 



L 2 diodo 7r V Lch 



\(f\H in t\i)\ 5 



(B13) 

(here bar denotes averaging with respect to polarization 
of both incident light and scattered quanta). Choose 
an initial ket-state |i) such that it contains an incident 
photon with momentum k;, frequency Qi = cki, and 
polarization Ai, and does not contain any QCB plasmon. 
This state can be written as |i) = |ki); |0) p , where 
|0) p is the plasmon vacuum, |kj)j = c^. A .|0);, and |0)j 

is the photon vacuum. A final bra-state (/| contains 
a scattered photon with momentum kf, frequency 
Qf = ckf, and polarization A/. It contains also a QCB 



plasmon P with the frequency uip (its quantum numbers 
will be specified below). The final state is written as 



(/| =(P| 0(^1, where (k y | = U f>X/ \0) 



A matrix element of the interaction which enters Eq. 
(EH, is 



where 



(f\H int \i) = (P|iJ int |0) p , 



Hint = (k/, Xf \H int \ ki, Ai) 



(B14) 



WgNL fev F \ 2 M-^ Qj 



P, 



k/ kj_ 

kf k{ 



(B15) 



In the case of normal incidence Eqs. I|B13|) - l|B15|) are 
equivalent to Eqs. © - l(T^|) . 
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